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Abstract. The phase structure of a simple Nambu-Jona-Lasinio model 

has been investigated at non-zero values of fi and H, where H is an external 

OS - 

Q"^ . magnetic field and /U is the chemical potential. On this basis magnetic oscilla- 

C^ . 

tions effects were considered. It was shown that there are standard (periodic) 

> . 

O . van Alphen-de Haas magnetic oscillations of some thermodynamical quanti- 

m 
(N 



0^ 



ties, including magnetization, pressure and particle density in the NJL system. 
Besides, we have found non-standard, i.e. non-periodic, magnetic oscillations. 



. since the frequency of oscillations is a ff-dependent quantity. Finally, there 

, arises an oscillating behaviour not only for thermodynamical quantities, but 



also for a dynamical quantity like the quark mass. 

Magnetic oscillations effects are well-known phenomena in condensed matter physics. In 
particular, the oscillation effect of the magnetization, which is called now the van Alphen-de 
Haas effect, was for the first time predicted by Landau and then experimentally observed in 
' some non-relativistic systems (in metals) more than sixty years ago |1|J^. At present, a lot 
of the attention of researchers dealing with magnetic oscillations is focused on relativistic 
condensed matter systems (mainly on QED at nonzero values of the chemical potential /i and 
external magnetic field H), since the results of these studies may be applied to cosmology, 
astrophysics and high energy physics 

It was shown in the framework of QED that the thermodynamical potential fi(yU, H) 
of the system has in 1-loop approximation the following form Q{fi,H) = Qmon{f^,H) + 
floscifJ', H), where fimon(/^5 H) is the monotonic part of f2(/x, if), and all magnetic oscillations 
are contained in the so-called oscillating part 
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nosc{fi,H) = J2[Ak{H) cos{2nku;) + Bk{H) sm{27rku;)], (1) 

k=l 

where u = (/i^ — m^)/(2eif) (e,m are electric charge and mass of fermions, respectively), 
and Ak{H), Bk{H) are smoothly varying functions. Due to the presence of trigonometric 
functions, expression (0) obviously oscillates over the variable {2eH)~^ with the frequency 
(/i^ — m^), which is not an if-dependent quantity. In condensed matter physics such kind 
of oscillations are usually called periodic ones. 

In the present talk magnetic oscillation effects are considered in the framework of quan- 
tum field theory with four-fermion interactions 

N ^ Q N N 

L = J2 Qkidqk + 77^[{J2 ^kqkf + {J2 Qki^bQk)^], (2) 

fc=l '^■^^ k=l k=l 

which is the A^-fermionic extension of the simplest Nambu - Jona-Lasinio model (NJL) ^ . [| 
Obviously, the model @ is invariant under (global) SU{N) and U{l)v transformations as 
well as continuous U{1)a chiral transformations: — ^ G^^"'''qk {k = 1, ...,N). 

We shall find the thermodynamic potential H), which is related to the corresponding 
effective potential Vff^t(S) of the NJL system (H) by 

ni^i,H) = V,Hi^)\ . (3) 

and contains all the information about thermodynamical quantities such as magnetization, 
particle density, etc. In the relation (^, one should first of all calculate the effective po- 
tential Vff^(S). So, before considering the magnetic oscillations, we can study the vacuum 
properties of the NJL model. 

Notice that special attention has been paid to the analysis of the vacuum structure of 
NJL-type models at non-zero temperature and chemical potential in the presence of 
external (chromo-)magnetic fields [3-[T0[|, with allowance for curvature and non-trivial space- 
time topology [pT] , p]2| . The combined influence of external electromagnetic and gravitational 



For simplicity, we consider in the following fermions ("quarks") of equal electric charge. 
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fields on the dynamical chiral symmetry breaking (DCSB) effect in four-fermion field theories 
was investigated in [0,|T^. However, the influence of both an external magnetic field H and 
chemical potential yU on the phase structure of the NJL model was not considered up to now. 

Phase structure of the model. The necessary information about the phase structure 
of a given field theoretical model is contained in the global minimum point of the corre- 
sponding effective potential. In the presence of /x, H the effective potential Vff^(S) of the 
NJL model has in leading order of large N the following form 



eH °° f 



A;=0 



Sk 



(4) 



where = 2 — 6ko, Sk = + 2eHk. VniX') is the effective potential at yU = 0, if 7^ 

Vh{^) = ^ + - ^{C'(-l, x) -\[x'-x]\nx + (5) 

where x = T?/{2eH), C,{v,x) is the generalized Riemann zeta-function, l,x) = 
dC,{u,x) I duly^^i, and 

is the effective potential at H, fi = 0. In (j^) A is the ultraviolet cut off parameter. Finally, 
let us remark that S is an auxiliary scalar field, which, at the tree level, is proportional to 
qq by the equations of motion. The global minimum point of the potential (|^) defines the 
vacuum expectation value of S and is equal to the dynamical quark mass. 

At fi, H = and G < Gc = Atc^/A"^ the global minimum point of Vo(S) equals to the 
value S = 0. Hence, in this case quarks are massless and chiral symmetry remains intact. 
If G > Gc, the effective potential (j^) has a nontrivial global minimum point, which we shall 
denote as M. (Evidently, M depends on the values of G and A [0].) 

At fi = 0, H ^ the chiral symmetry of the model is spontaneously broken for arbitrary 
values of the bare coupling constant G. This is due to the fact, that the global minimum 



point T,o{H) of the potential Vh(E) is unequal to zero [§|,|10 . 

In order to study the properties of the NJL model vacuum for the general case, when 
both fi and H are nonzero, one should find all solutions of the stationarity equation 



Sk 







(7) 



and select that one, at which the potential Vh'^(S) takes its smallest value. This is the global 
minimum point for the function (^. The properties of this point as a function of fi and H 
give us a lot of information about the ground state. We omit here the detailed consideration 
of this procedure and present directly the phase structure description of the model (Figure 
1). 




FIGURE 1. Phase structure of the NJL modeL (Detailed description of the figure is given in the text.) 



In this figure, in the plane (/i, y/eH) the phase portrait of the model is qualitatively 
represented for the case Gc < G < (1.225...)Gc, where M is the quark mass at fj, = H = 0, 
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Ml = (A^/2 — 27r^/G')-^/^. Here one can see infinite sets of symmetric massless Aq,Ai,... 
pliases, as well as massive phases Cq, Ci, ... with DCSB. In addition, there is another massive 
phase B. Dashed and solid lines in Figure 1 are critical curves of first- and second-order 
phase transitions, respectively. One can also see on this phase portrait infinitely many 
tricritical points tfc,Sfc {k = 0,1,2,...) which lie on the boundary between massless and 
massive phases (chiral boundary). (A point of the phase diagram is called a tricritical 
one if, in an arbitrarily small vicinity of it, there are first- as well as second-order phase 
transitions.) Numerical investigation gives the following values of the external magnetic 
field corresponding to tricritical points to and sq at different values of the bare coupling 
constant G: eHtJA^=0.01...; 0.08...; 0.13... as well as eHsJA^=0.006...; 0.056...; 0.103... 
for G/Gc=l.01; 1.1; 1.2, respectively. We should also remark that the part tondH) of the 
chiral boundary is described by the equation Vff^(O) = Vff^(So (//)). 

Points {fi, H) of the phase diagram, lying above the chiral boundary, correspond to the 
chirally symmetric ground state of the NJL model. One-fermion excitations of this vacuum 
have zero masses. At first sight, it might seem that the properties of this symmetric vacuum 
are slightly varied, when parameters /i and H are changed. However, this is not the case, 
and in this region, as was mentioned above, we have infinitely many massless symmetric 
phases of the theory corresponding to infinitely many Landau levels, as well as a variety of 
critical curves of second-order phase transitions. Let us next show this. 

It is well-known that the state of thermodynamic equilibrium (the ground state) of an 
arbitrary quantum system is described by the thermodynamic potential (TDP) Q, which is 
just the value of the effective potential at its global minimum point (see (§)). In the case 
under consideration, the TDP Q{fj,, if ) at /x > Mi (see Figure 1) has the form 



where = y2eHk. We shall use the following criterion of phase transitions: if at least one 
first (second) partial derivative of fi(yU, H) is a discontinuous function at some point, then 
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this is a point of a first- (second-) order phase transition. 

Using this criterion, let us show that hnes /fc = {(/i, if ) : = \/2eHk} {k = 1,2,...), are 
critical lines of second-order phase transitions. 

Indeed, from (^) one easily finds 



an 



on 



(9) 



as well as: 



eHfj, 



-oo. (10) 



Equation (|^) means that the first derivative dQ/dfi is a continuous function on all lines Ik- 
However, the second derivative d'^Q/{dfiY has an infinite jump on each line Ik (see (p!OD), so 
these lines are critical curves of second-order phase transitions. (Similarly, we can prove the 
discontinuity of d^il/{dHY and d'^Vt/dfJidH on all lines /„.) 

The presence of an infinite set of massive phases Ck on the phase portrait is conditioned 
by a special structure of the stationarity equation (|^. Analytical and numerical considera- 
tions of it show that below the chiral boundary the effective potential global minimum point 
S(yU, ii), which is identical to the quark mass, has /i and H dependences. The function 
S(/i, if) is a continuous one inside each of regions Ck- However, it is a discontinuous one 
on each of the curves Mt^, where the quark mass changes its value by a jump. That is 
why boundaries between C^-regions are the first order phase transition lines. In contrast, 
in the phase B, the global minimum point is equal to T,q{H) (= quark mass in the case 
/i = 0,ii 7^ 0), which is a //-independent quantity. This means that the particle density 
n = —dQ/d^ in the ground state of the phase B is identically equal to zero, whereas in each 
phase Ck this quantity differs from zero. 

Magnetic oscillations. Now we want to show that there arise, from the presence of 
infinite sets of massless Ak phases as well as of massive Ck ones, magnetic oscillations (the 
so-called van Alphen-de Haas-type effect) of some physical parameters in the NJL model 
gauged by an external magnetic field. 
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Let the chemical potential be fixed, i.e. fi = const > Mi (see Figure 1). Then on the 
plane (/i, VeH) (see Figure 1) we have a line that crosses the critical lines h, ... at points 
corresponding to some values Hi,H2, ... of the external magnetic field. The particle density 
n and the magnetization m of any thermodynamic system are defined by the TDP in the 
following way: n = —dQ/dfi, m = —dQ/dH. At fi = const these quantities are continuous 
functions of the external magnetic field only, i.e. n = n{H), m = m{H). We know that all 
the second derivatives of fl{fi,H) are discontinuous on every critical line /„ (see (0)). The 
functions n{H) and m{H), being continuous in the interval H G (0, oo), therefore have first 
derivatives that are discontinuous on an infinite set of points Hi, H^, ... Such a behaviour 
manifests itself a phenomenon usually called oscillations. 

Analogously to QED and condensed- matter physics let us again separate the ex- 

pression for a physical quantity with oscillations into two parts: the first monotonic one 
does not contain any oscillations, whereas the second part, which is of particular interest 
here, contains all the oscillations. Following this rule, we can write down, say, the TDP (^ 
of the NJL model in the form Q{fi,H) = Qmonif^, H) + QoscifJ', H). In order to present the 
oscillating part QoscifJ, H) in an analytical form, we shall use the technique elaborated in [Q, 
where manifestly analytical expressions for this quantity were found in the case of a perfect 
relativistic electron-positron gas. This technique can be used without any difficulties in our 
case, too. So, applying in (§) the Poisson summation formula [|l| 

oo oo °? 

^Q;„$(n) =2^afc / $(x) cos(27r/cx)c?a;, (11) 

n=0 fc=0 

where a„ = 2 — 5„o; one can get for Vtosci.lJ'i H) the following expression 

^""^ = 51 [Qiirkiy] cos(7rA;i/ + n/A) + P{TTkiy) cos(7rfcz/ - 7r/4)], (12) 

where u = fi^/{eH). The functions P{x) and Q{x) in (|T2|) are connected with the Fres- 
nel integrals C{x) and S{x) [|T6|: C{x) = | + \f^[P{^) ^^^^ + Q{x) cosx], S{x) = 
' ^[P{x) cosx — Q{x) sinx]. They have, at a; — > cxd, the following asymptotics |]16 



P{x) = X ^—3x ^ /A + Q{x) = —X ^/2 + 15x ^/8 + ... Formula ([T^) presents, in a mani- 
festly analytical form, the oscillating part of the TDP (||) for the NJL model at /i > Mi. In 



the case under consideration, since the TDP is proportional to the pressure of the system, 
one can conclude that the pressure in the NJL model oscillates when H —>■ 0, too. It follows 
from (|1^) that the frequency of oscillations over the parameter {eH)~^ equals and does 
not depend on H. So, in this case we have periodic magnetic oscillations. Then, starting 
from (|T^), one can easily find the corresponding expressions for the oscillating parts of n{H) 
and m{H). These quantities oscillate at if ^ with the same frequency and have a 
rather involved form, so we do not present them here. 

Finally, we should note that the character of magnetic oscillations in the NJL model 
at yU > Ml resembles the magnetic oscillations in massless quantum electrodynamics [^Q. 
Indeed, in this case in both models one can find periodic magnetic oscillations of some 
thermodynamic parameters. 

Now let us show that at a fixed value of the chemical potential and M < /i < Mi the 
character of magnetic oscillations is changed. In this case on the plane (/i, y/eH) we have 
a line drawn through an infinite set of the C^-phases. Hence, the thermodynamic potential 
of the NJL system has the following form: Q{fi,H) =Ve^^(S(yU, if)), where ii) is the 
global minimum point of the potential Vi^^(S). Applying in (4) again the formula (0), one 
can find the following expression for the oscillating part of TDP 



°° 'el 



3/2 

f^osc-y^f^ [Q(7rA;z/)cos(27rA;c^ + 7r/4) + P(7rA;z/)cos(27rA;u; -7r/4)], (13) 



k=l 

where z/ = fi'^/{eH), u = {fi'^ - E'^{fi, H))/{2eH). From one can see that the TDP 
H) oscillates with frequency (yU^ — H))/2 if the variable (eH)^^ tends to infinity. 
Since H) is, up to a sign, equal to the pressure in the ground state of the system, 
also in the present case the pressure in the NJL model is an oscillating quantity. Moreover, 
other thermodynamic quantities such as particle density n = —dQ/dfi and magnetization 
m = —dVt/dH oscillate with the same frequency. 

Here we should do an important remark. In the NJL model at M < < Mi, in 
contrast to QED, the magnetic oscillation frequency is a ii-dependent quantity. (Since the 
quark mass T,{fi,H) has ii-dependency. ) So, strictly speaking, in the NJL model magnetic 
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oscillations are not periodic ones. Recently, similar peculiarities of magnetic oscillations 



are observed in some ferromagnetic semiconductive materials such as HgCr2Se4 [jT^, where 
non-periodic magnetic oscillations over the variable {eH)^^ were found to exist for electric 
conductivity as well as magnetization. 

Finally, we should remark that in the NJL model not only thermodynamic quantities 
oscillate, but some dynamical parameters of the system do as well. This concerns, in partic- 
ular, oscillations of the dynamical quark mass. In fact, by applying the Poisson summation 
formula (0) to the stationarity equation (0) and searching for the solution S(/i, if) of this 
equation in the form S(yU, if) = S^on + ^osc, one can easily find the following expressions 
for H ^0: 

{eHf/^ ^ sin(27rA:^ - 7r/4) 

^OS.{,, H) r. g — , (14) 

where a) = (/i^ — M^) / {2eH), and M = M(/i) is the quark mass at H = 0, fx ^ 0. 

Conclusions: Let us point out once more that for strongly correlated fermionic systems 
there is a possibility to observe nonperiodic magnetic oscillations. Moreover, in such systems 
in the presence of an external magnetic field some dynamical quantities (for example, fermion 
masses) should oscillate as well. Our results may be applicable in astrophysics, in the physics 
of neutron stars etc, where one should take into account the relativistic character of different 
phenomena. 

Note that the strength of the surface magnetic field of a neutron star is about 10^^ G 



and in the interior it is probably 10 G |18|. Our numerical estimates of the ii^g values 



using A = 700 Mev show that the magnetic field corresponding to the tricritical point sq 
varies in the interval 10^^ G^ 10^^ G, when 1.01 < G/Gc < 1.2. Hence, the typical neutron 
star magnetic field strengths are much smaller, than the value of Hg^, and are located in 
the oscillation region of the NJL model (see Figure 1). So, the if-dependency of different 
physical parameters (such as particle density, magnetization, quark mass, etc) inside neutron 
stars possibly has a nonperiodic oscillating character. 

Despite the relativistic character of our investigations, we believe that qualitatively the 



presented results are valid for nonrelativistic electronic systems, and may be applicable in 
condensed matter physics, too. 

More complete information about phase structure as well as magnetic oscillations in 



several NJL-type models one can find in our recent paper ||T9| 
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